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We calculate the 3-loop bare /3-function of QCD, formulated on the lattice with the clover fermionic action. The 
dependence of our result on the number of colors N, the number of fermionic flavors Nf, and the clover parameter 
csw, is shown explicitly. A direct outcome of our calculation is the two- loop relation between the bare coupling 
constant go and the one renormalized in the MS scheme. Further, we can immediately derive the three-loop 
correction to the relation between the lattice A-parameter and go, which turns out to be very pronounced. 

We also calculate the critical value of the hopping parameter, acc, in the clover action, to two loops in pertur- 
bation theory. This quantity is an additive renormalization; as such, it exhibits a linear divergence in the lattice 
spacing. We compare our results to non perturbative evaluations of Kc coming from MC simulations. 



1. INTRODUCTION 

The clover action for lattice fermions was in- 
troduced a number of years ago |^ , as a means of 
reducing finite lattice spacing effects, ft is widely 
used nowadays in Monte Carlo simulations. 

To monitor the onset of the continuum limit, 
tests of scaling must be performed on measured 
quantities. In particular, asymptotic scaling is 
governed by the bare /^-function: 

/?L(ffo) = (1) 

(a is the lattice spacing, g (go) the renormalized 
(bare) coupling constant, fj, the renormalization 
scale). For go one may write (3l as: 

/5l(5o) --feo.9o-^iffo-^25o + -- (2) 

The first two coefficients, 6o and 5i, are univer- 
sal and well known in SU{N) gauge theory with 
Nf fermion species; the 3-loop coefficient 63 , on 
the other hand, is regularization dependent. In 
the case at hand, 62 is thus expected to depend 
not only on N and Nf , but also on the parameter 
Csw of the clover action (see next Section). 
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We calculate fof for arbitrary N, Nf and csw- 
The analogous calculation for pure gauge the- 
ory without fermions H,^, as well as for Wilson 
fermions Q], was done a few years ago. We follow 
the general setup of those publications. 

The /3-function enters directly into the relation 
defining the parameter A^, : 

■ [1 + 995 + ...], 5 = (if - iiodf ) / 24=. (3) 

The "correction" factor q turns out to be very 
pronounced for typical values of csw and 50- 

A direct outcome of our calculation is the 
two-loop relation between the MS coupling a = 
5V(47r) and ao = (?g/(47r) : 

a = ao + di(a/i) (^0 + ^2(0^1) al + O (ag) , (4) 

This relation is useful in studies involving running 
couplings or renormalized quark masses. 

In Sec. ^ we present our results for 63 , di(a/i) 
and d2(a/Lt), as functions of N, Nf and csw- Fur- 
ther technical details and checks of our calcula- 
tions are relegated to a longer write-up 

We have also calculated the critical value of the 
hopping parameter Kc, to two loops, using the 
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clover action. Since Wilson fermions break chi- 
ral invariance explicitly, merely setting their bare 
mass to zero does not ensure chiral symmetry in 
the continuum limit; quantum corrections intro- 
duce an additive renormalization to the fermionic 
mass, which must then be fine tuned to a vanish- 
ing renormalized value. Consequently, the hop- 
ping parameter k is shifted from its naive value. 

The additive mass renormalization is linearly 
divergent with the lattice spacing. This adverse 
feature of Wilson fermions poses an additional 
problem to a perturbative treatment. Indeed, our 
calculation serves as a check on the limits of ap- 
plicability of perturbation theory, by comparison 
with non perturbative Monte Carlo results. 

In the present work we follow the procedure of 
Rcf . 1^ , in which Hc was computed using Wilson 
fermions without 0{a) improvement. In Sec. |^ we 
present our results on Kc, showing explicitly the 
dependence on N, Nf and csw- Details on our 
calculation can be found in our publication [Q . 

2. THE P FUNCTION 

Our starting point is the Wilson formulation of 
the QCD action on the lattice, with the addition 
of the clover Q fermion term, Ssw: which reads 
in standard notation: 

5 _ 



(5) 



S'sw — Csw ^ 

X, fJ., iJ, f 

[Ffj.i' = g {Q ^ Qufi), 

Here C/^^ ^ (x) is the usual product of link variables 
(x) along the perimeter of a plaquette in the fi- 
V directions, originating at a;; / is a flavor index; 
f/ji/ = (z/2)[7^, The value of the parame- 
ter Csw can be chosen arbitrarily; it is normally 
tuned in a way as to minimize 0(a) effects. The 
lattice /^-function is independent of the renormal- 
ized fermionic masses, which may be set to zero. 

We compute the relation between go and de- 
fined in the MS renormalization scheme: 



50 = ^g(3o, a[i)g. 



(6) 



The one- and two-loop terms of have the form: 
Zg(5o, a,[if = 1 + LQ{a^i)gl + Li{a^) + ... 



Lo{x) ^ 2bo\nx + Iq, Li{x) ^ 2bi\nx + h. (7) 

The constant Iq is related to the ratio of the 
associated A parameters: 

?o = 26oln(AL/Aj^) . (8) 

Its value is known (see e.g. Ref. Q and references 
therein) and is presented here with increased ac- 
curacy for the csw"dependent coefficients: 

Iq^ 0.1699559991998031(2) TV + iV/ /qi 

loi = 0.006696001(5) - csw 0.00504671402(1) 
-fcsw^ 0.02984346720(1) (9) 

The dependence on csw is quite pronounced, 
leading to changes in A/, of up to a factor of 2. 



The quantity 63 can be obtained from Iq, I 
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6^ = 62 



bik + boh. (10) 

where &2 is known from the continuum. Com- 
puting li amounts to a two- loop calculation of 
the one-particle irreducible two-point function of 
a background gauge field. We find: 



hi = 



I12 = 



-3/(128Af^) -t- Nf [ hi /N + h2 N] 
-1-0.018127763034(4) - A^^ 0.0079101185(2) 
-0.0011967(14) + csw 0.0001578(3) 
+csw^ 0.0052931(2) + csw^ 0.00050624(3) 
-Hcsw'' 0.00008199(1), 
0.0009998(16) + csw 0.0000342(4) 
-csw^ 0.0048660(6) - csw^ 0.00021431(3) 
-csw"^ 0.00004382(1). (11) 

Substitution in Eq. (10) now yields &2 > 
value of N , Nf, csw- The correction coefficient 
q of Eq. (H) also follows immediately; it brings 
about a substantial correction to asymptotic scal- 
ing, with a pronounced csw dependence. Finally, 
the coefficients di {a/i) and d2 (a/i) can be read off 
Eqs. (0), (i), (I), (0). 

There exist several constraints on the alge- 
braic and numerical values of individual dia- 
grams. A particularly strong check is provided 
by Ref. Eq. (5.6) in that reference reads: 
d2(A^=3, c^°^=l) = 1.685(9) - 8.6286(2) c^^^. 



For the same quantity, our results lead to: 
1.6828(8) -8.62843775(1) c^\^. Both sets of num- 
bers are clearly in very good agreement. 
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3. THE HOPPING PARAMETER 

We set out to calculate the hopping parameter, 

K = 1/ (2TOBa + 8r), (12) 

which is an adjustable quantity in numerical sim- 
ulations; TOs is the bare fermionic mass, and r is 
the Wilson parameter appearing in the fermionic 
action, usually set to 1. The critical value of k, 
at which chiral symmetry is restored, is thus l/8r 
classically, but gets shifted by quantum effects. 
For a vanishing renormalized mass we require: 

TOB = E-^(0,m_B,.go) (13) 

where I]^(p, tob, go) is the truncated, IPI 
fermionic two-point function. We solve this recur- 
sive equation for uib perturbatively. We write: 

5]^(0, mB,go) = gl + ^^^'^ + ' ' ' (14) 
Two diagrams contribute to S'-'^', and 26 to 
S*^^). Certain sets of diagrams must be evalu- 
ated together for infrared convergence. The de- 
pendence on cgw is polynomial. For E^^^ we find: 

= {N^ -I)/ N (-0.16285705871085(1) (15) 

+CSW 0.04348303388205(10) 

+csw^ 0.01809576878142(1) ) 

One- and two-loop results pertaining to csw = 
are as in Ref. and can be found with greater 
accuracy in a subsequent work |^ . 

For csw 7^ 0, only one-loop results exist so 
far in the literature; a recent presentation for 
Ar=3, csw=l @ agrees with our Eq. (|l5|): 

^ -0.2700753495(2), Ref. 

-0.2700753494597(5), Eq. (|l5|). 

We now turn to the much more cumbersome 
evaluation of E'^'. Our result is: 

E(2)/(7V2-1) = 

[(-0.017537(3) -I- 0.016567(2) 

+Nf/N 0.00118618(8)) 
+ (0.002601(2) - 0.0005597(7) 

-Nf/N 0.0005459(2)) csw 
-H(-0.0001556(3)-h 1/A^2 0.0026226(2) 

+Nf/N 0.0013652(1)) csw^ 
+ (-0.00016315(6)+ 1/A^2 0.00015803(6) 

-Nf/N 0.00069225(3)) csw^ 
+ (-0.000017219(2)+ l/iV2 0.000042829(3) 

-Nf/N 0.000198100(7)) csw''] 



(16) 



Several important consistency checks can be 
performed on the values of individual diagrams; 
our results satisfy these checks . Eqs. (|l^, [l^ ) 
lead immediately to the critical mass and hopping 
parameter, making use of Eqs. (|3|, |l2|) . 

A number of non-perturbative determinations 
of Kc exist for particular values of go, csw, and 
Nf . We present these in Q , along with our results 
and with "dressed" results obtained by a resum- 
mation of cactus diagrams |l^ . Comparing with 
the Monte Carlo estimates, dressed results show 
a definite improvement over non-dressed values. 
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